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ABSTRACT 

We show how the continuity equation can be used to determine pattern speeds in the 
Milky Way Galaxy (MWG). This method, first discussed by Tremaine & Weinberg 
in the context of external galaxies, requires projected positions, (l,b), and line-of- 
sight velocities for a spatially complete sample of relaxed tracers. If the local standard 
of rest (LSR) has a zero velocity in the radial direction (ulsr)j then the quantity 
that is measured is AV = fl p Ro — Vlsr, where tt p is the pattern speed of the non- 
axisymmetric feature, Ro is the distance of the Sun from the Galactic centre and 
Vlsr is the tangential motion of the LSR, including the circular velocity. We use 
simple models to assess the reliability of the method for measuring a single, constant 
pattern speed of either a bar or spiral in the inner MWG. We then apply the method 
to the OH/IR stars in the ATCA/VLA OH 1612 MHz survey of Sevenster et al., 
finding AV = 252 ± 41 km s -1 , if ulsr = 0. Assuming further that Rq — 8 kpc and 
Vlsr = 220 km s -1 , this gives fl p = 59 ± 5 km s -1 kpc -1 with a possible systematic 
error of perhaps 10 km s -1 kpc -1 . The non-axisymmetric feature for which we measure 
this pattern speed must be in the disc of the MWG. 

Key words: Galaxy: bulge - Galaxy: kinematics and dynamics - Galaxy: disc - 
Galaxy: structure - galaxies: spiral 



1 INTRODUCTION 

Our position within the Milky Way Galaxy (MWG) has 
made study of its large-scale structure and dynamics dif- 
ficult. This is partly due to the Sun's position within the 
dust layer of the MWG, as well as the difficulty of deter- 
mining distances. Nevertheless, in recent years, considerable 
progress has been made in understanding our Galaxy. 

It is now clear that the MWG is a barred galaxy, as was 
first proposed by de Vaucouleurs (1964). Evidence of this 
comes from the COBE/DIRBE near infrared light distribu- 
tion (Weiland et al. 1994; Dwek et al. 1995; Binney et al. 
1997), star count asymmetries (Nakada et al. 1991; White- 
lock & Catchpole 1992; Sevenster 1996; Nikolaev & Wein- 
berg 1997; Stanek et al. 1997; Hammersley et al. 2000), gas 
dynamics (Peters 1975; Cohen & Few 1976; Liszt & Burton 
1980; Gerhard & Vietri 1986; Mulder & Liem 1986; Binney 
et al. 1991) and the large microlensing optical depth towards 
the bulge (Paczynski et al. 1994; Zhao et al. 1996). A recent 
review of the structure of the bulge and disc can be found 
in Gerhard (2001). 

The principle dynamical parameter of a bar is its pat- 
tern speed, fl p , since it determines the orbital structure of 
stars in the disc. Bar pattern speeds can be parametrized by 
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the distance- independent ratio 7Z = Dl /olb, where Dl is the 
Lagrangian radius at which the gravitational and centrifu- 
gal forces cancel out in the bar rest-frame (Dl is usually 
approximated by the axisymmetric corotation radius) and 
as is the bar semi-major axis. A bar is termed fast when 
1.0 < 1Z < 1.4. Knowledge of bar pattern speeds constrains 
the dark matter (DM) content of disc galaxies: Debattista 
& Sellwood (1998, 2000) argue that fast bars require that 
the disc they inhabit must be maximal, in the sense that the 
luminous disc provides most of the rotational support in the 
inner galaxy. Some evidence for fast bars comes from hy- 
drodynamical models of gas flow, particularly at the shocks. 
Three such studies are: NGC 1365 (7Z = 1.3, Lindblad et al. 
1996), NGC 1300 (71 = 1.3, Lindblad & Kristen 1996), and 
NGC 4123 (71 = 1.2, Weiner et al. 2001). 

In the MWG, Q p can be estimated by comparing hy- 
drodynamical models of the inner MWG with the observed 
longitude- velocity, (l,v), diagram. Several such studies have 
been carried out: Englmaier & Gerhard (1999) found Dl = 
3.5 ± 0.5 kpc, Q p = 59 ± 2 km s _1 kpc" 1 and 71 = 1.2 ± 0.2, 
while Fux (1999) found D L = 4 - 4.5 kpc and Q p = 35 - 45 
km s" 1 kpc -1 and the preferred model of Weiner & Sell- 
wood (1999) had D L = 5.0 kpc, Q p = 41.9 km s _1 kpc" 1 
and 7Z = 1.4. Binney et al. (1991) interpreted the (l,v) dia- 
gram for HI, CO and CS emission using orbits in a barred 
potential with Q p — 63 km s _1 kpc -1 . Another method 
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for measuring Q p in the MWG involves the identification 
of the action of resonances: Binney et al. (1997) suggested 
that a local density maximum along the minor axis seen 
in the dust-corrected infrared COBE/T)IRBE data could be 
identified with the location of a Lagrange point, which gave 
Sip = 65 ± 5 km s _1 kpc -1 , having assumed TZ ~ 1.2 and 
given the value of as from the deprojection. Dehnen (1999) 
argued that the bimodal distribution of stellar velocities 
in the solar neighborhood, as observed by HIPPARCOS, is 
due to the action of the outer Lindblad resonance (OLR), 
whereby he obtained fi p = 53 ± 3 km s -1 kpc -1 . Sevenster 
(1999) suggested that features in the distribution of OH/IR 
stars can be used to locate the corotation and inner Lind- 
blad resonances, giving Q p ~ 60 km s _1 kpc -1 (TZ <; 1.4). 
The differences in these pattern speed measurements reflect 
uncertainties in viewing geometry, as, rotation curve and 
bar axis ratio. 

The MWG also has spiral structure in it, which is less 
well constrained. Even the number of arms, whether two 
or four, continues to be an issue of some discussion. Vallee 
(1995) reviews a variety of different observational tracers, 
and concludes that the spirals are four-armed, logarithmic 
and have a pitch angle, p ~ 12° ± 1°. More recently, Drim- 
mel & Spergel (2001) found that a two arm spiral structure 
dominates the non-axisymmetric near-infrared emission. 

The pattern speeds of spirals are generally quite poorly 
constrained. When a bar is present, as in the MWG, the 
expectation that spirals are driven by it (e.g. Sanders & 
Huntley 1976) would seem to require that the bar and spi- 
rals share a common pattern speed. This idea appears to be 
further reinforced by the observation that many such spi- 
rals connect with the bar ends. Sellwood & Sparke (1988), 
however, presented iV-body simulations in which bars and 
spirals had different pattern speeds, but with the spirals still 
generally connected to the bar ends. Tagger et al. (1987) and 
Sygnet et al. (1988), proposed that the pattern speeds of the 
bar and spirals are non-linearly resonantly coupled, with the 
corotation of the bar marking the inner Lindblad resonance 
of the spiral. Toomre's (1981) swing amplification theory of 
spiral structure formation requires a radially (and possibly 
temporally) varying pattern speed. Modal theories of spiral 
formation, on the other hand, depend on a (nearly) constant 
spiral pattern speed in space (but not necessarily in time, see 
Bertin & Lin 1996). Measurements of spiral pattern speeds 
in external galaxies are still quite few, generally assume a 
radially constant £l p and tend to give conflicting results. For 
example, the spiral pattern speed in M81 has been measured 
in a variety of ways but so far no unique pattern speed has 
been found (see discussion in Canzian 1993). 

The situation is not much better in the MWG. Mor- 
gan (1990) was able to model the distribution of pulsars 
within roughly 10 kpc of the Sun with spirals of fl p — 13.5 
km s _1 kpc -1 , regardless of the arm multiplicity. On the 
other hand, Amaral & Lepine (1997), using the positions 
and ages of young open clusters within ~ 5 kpc of the Sun, 
estimated fi p ~ 20 km s -1 kpc -1 . Mishurov & Zenina (1999) 
used the velocity field of a sample of Cepheids which were 
less than 5 kpc from the Sun to conclude that the Sun is 
displaced outwards from the corotation by ~ 0.1 kpc (then 
£l p ~ J1 Q ~ 27 km s" 1 kpc" 1 ). 

A direct method for measuring S7 P in systems satisfying 
the continuity equation was developed by Tremaine & Wein- 



berg (1984). This method has successfully been applied to 
four barred galaxies so far: NGC 936 (TZ — 1.4 ±0.3, Merri- 
field & Kuijken 1995), NGC 4596 (TZ = 1.15±g"J|, Gerssen et 
al. 1998), NGC 7079 (TZ = 0.9 ±0.15, Debattista & Williams 
2001) and NGC 1023 (TZ = 0.77±°;22, Debattista et al. 2002). 
In view of the spread in previous measurements of pattern 
speeds in the MWG, it is worth considering whether this 
method may not also be applied it. Kuijken & Tremaine 
(1991) have already derived an equivalent of the Tremaine- 
Weinberg method for the MWG; their main goal was to con- 
strain the radial motion of the local standard of rest (LSR) . 
Applying this method to the HI distribution, they found 
that pattern speeds from Q, p — to fl@ are possible for an 
LSR with a radial velocity, mlsr = —4 km s _1 to 8 km s . 

In this paper, we explore the Tremaine- Weinberg 
method for the MWG. This method, unlike others, has the 
considerable advantage of being model independent. In the 
general case, the only assumption made is that the non- 
axisymmetric density distribution is in stationary cylindri- 
cal rotation (i.e. Sl p = S} p (p), where p is the distance from 
the Galactic centre). While this assumption breaks down 
in some cases (e.g. if the MWG disc is strongly interacting 
with a companion) , it is expected to be a very good approx- 
imation inside the solar circle. Since we also wish to apply 
the method to real data, we have made a number of sim- 
plifying assumptions. For example, we assume that there is 
one unique pattern speed in the inner MWG; if this is not 
the case, then some density and asymmetry weighted aver- 
age pattern speed is measured. The value of Q p obtained 
depends sensitively on the assumed radial velocity of the 
LSR; however this dependence can be stated explicitly. The 
method requires complete samples of tracers: one such sam- 
ple is obtained from the ATCA/VLA OH 1612 MHz survey 
of OH/IR stars. 

In we derive the three-dimensional version of the 
Tremaine- Weinberg method for the MWG, with the solar 
and LSR motion taken into account. We test the method in 
using simple models, exploring under what circumstances 
the method works best. In we introduce the ATCA/VLA 
OH 1612 MHz survey, and describe how we selected a sam- 
ple of older, relaxed OH/IR stars from it. We also discuss 
our completeness corrections for the survey data. Then in 
S^, we present the pattern speed analysis for our sample of 
OH/IR stars. Since the integrals in the Tremaine- Weinberg 
method measure asymmetries between positive and negative 
longitude, they are susceptible to noise and/or observational 
bias; therefore in we present corroborating evidence for 
the signal we find. Finally, in we discuss our results in 
terms of MWG structure and prospects for using the same 
method with future data sets. 



2 THE TREMAINE- WEINBERG METHOD 
FOR THE MILKY WAY 

Starting from the continuity equation and the assumption 
that the visible surface density, p, of a barred (or otherwise 
non-axisymmetric) disc depends on time, t, in the following 
simple way: 

v(x,y,t) = (i(p,ip -n. p t), (1) 
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where (x, y) are arbitrary cartesian coordinates in the plane 
of the disc, and (p, ip) the corresponding polar coordinates, 
Tremaine & Weinberg (1984) showed that the pattern speed, 
Q„, can be written as: 



/I! fe ( y )/-°°oc Vios(X,Y)n(X,Y)dXdY 
A 00.TL Xfi(X,Y)dXdY 



(2) 



Here, Vj os is the line-of-sight (minus the systemic) veloc- 
ity, h(Y) is an arbitrary continuous weight function, i is 
the inclination angle and (X, Y) are galaxy-centered coordi- 
nates along the disc's apparent major and minor axis respec- 
tively. The Tremaine- Weinberg (TW) method needs a tracer 
of the non-axisymmetric feature which satisfies the continu- 
ity equation. The old stellar population in undisturbed SBO 
type disc galaxies provides an ideal tracer of the surface den- 
sity. The short lifetimes of massive stars precludes their use 
as tracers. Moreover, gas cannot be used as a tracer because 
of its conversion between atomic and molecular states, al- 
though modelling may be used to describe such conversions 
(Bienayme et at. 1985). 



2.1 Two dimensions 

We now seek to derive an equivalent to the TW method for 
the MWG. Our position within the MWG disc gives a unique 
viewing geometry which does not permit the standard TW 
method of Equation (^) to be used. As shown by Kuijken & 
Tremaine (1991), however, it is still possible to use the con- 
tinuity equation to derive an expression for a pattern speed 
in the MWG. Here, we re-derive their two-dimensional ex- 
pression before considering the three-dimensional case. As- 
suming Equation (jj]) holds, the continuity equation in the 
plane of the disc can be written: 



dp dp 



Ox 



dy 



V ■ (pv) = 0. 



(•3) 



Switching to Sun-centered cylindrical coordinates (R,l), 
(where I is the usual Galactic longitude and R is the dis- 
tance from the Sun), this becomes: 



- - a 



dp f R cosl\ dp 

- Rosml dR + { 1 -—R—)m 
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1 d(Rpv f 



+ 



1 d(pvi) 



0. 



(4) 



R dR R dl 

In this equation, (vr,vi) is the velocity in the (R, I) frame 
and Ro is the distance of the Sun from the Galactic centre. 
Typically, the component vi is not available; fortunately, this 
term drops out after an integration over I: 



p sin Idl — 
P R 



p sin Idl 



+ 



RdR 



pvR.Rdl = 0. 



(5) 



If we now multiply by Rg(R), where g(R) is an arbitrary 
function vanishing at the Sun, and integrate over R, we ob- 
tain, after some rearrangement: 



QpRo 



JT JT vf( R >R dl dR 

JT.C 'M/CR)sinZ dl dR 
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where we have assumed that pg(R) goes to zero at infinity 
and f(R) = Rdg/dR. All the quantities on the right hand 
side of the above equation are observable. 

In the case of discrete tracers, the surface density gets 
replaced by a sum over delta functions and the above ex- 
pression is replaced by: 



QpRo 



J^JiR^smh 



(7) 



Thus f(R) can be thought of as a detection probability for 
the discrete tracers. 

In Equation (^) above, we integrated I from to 2ir. 
But if the non-axisymmetric part of p is limited to the range 
[— Iq, lo], it is then trivial to show that these integration lim- 
its are sufficient, since axisymmetric components do not con- 
tribute to the integrals, and all surface terms resulting from 
the I integrals vanish by symmetry. It is necessary, however, 
that the integrals extend from — lo to Iq (i.e. integration in 
the range [—li,h], with l\ 7^ I2 gives incorrect results.) 

When the bar is either along or perpendicular to the 
line of sight, the numerator and denominator of Equation 
([j]) are zero, and Q p Ro cannot be determined. This is to be 
compared with the standard TW equation, where fi p cannot 
be measured for a bar along the apparent major or minor 
axis of the disc. 



2.2 Three dimensions 

In 3D, the derivation is more or less the same. The Sun- 
centered spherical polar coordinates are now (r,l,b), with 
(l,b) being the standard galactic coordinates and r the dis- 
tance from the Sun. Following a triple integral much in the 
spirit of the 2D case, we obtain: 

f„°° f + , b ° f + . f(r) p cos b v r dl db dr 

Q ft _ JO J -bp J-lg ■> K I ^ ^ . . 

P ° JT C / W " cos2 b sin ldldbdr ' 

where 60 must be such that all the non-axisymmetric parts 
are included. If the system is symmetric about the mid- 
plane, then the b integral can extend over the region [0, ±60] 
only: data of this kind are easier to obtain. A discrete version 
of this equation is given by: 



QpRo 



i _ /C(Zo) 

J2i fi r i) sink cos bi V(lo) 



(9) 



where we have introduced a new notation for the integrals. 
(Note that one factor of cos b drops out because of the Jaco- 
bian in the delta functions^ In the limit that all the tracers 
are in a plane, Equation reduces to Equation (Q) of the 
2D case. 

Both V and K, are defined as integrals/sums over I rang- 
ing from — lo to Iq. Since sin I and v r change sign across / = 0, 
the integrals will partially cancel out. Thus V and /C largely 
measure a difference across I = 0. When the difference is 
small, V and /C are small and the measured fl p Ro will be 
poorly determined. This can be seen from the expression 



for the variance in f2 



+ 



To be more explicit by what we mean by small, we define 
V± = J °° J ±l ° f(r)pcos 2 bsinl dl db dr, with a simi- 
lar definition for K.±. (In this notation V = P+ — P_ and 
similarly for K. ) Then we introduce the variables 
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Air = 



|/C_| + ]^ 



A-p = 



|P| 



IP- 



(10) 



These two parameters, which take values in the range [0, 1], 
quantify the degree of global asymmetry in the system. 
When A-p — A>c = 0, no asymmetry signal is present (but 
note that the system need not be axisymmetric, or even sym- 
metric in projection) while when A-p = A/c = 1 the observed 
system is maximally asymmetric. In ^ we use simple mod- 
els to explore the dependence of measurement errors on A-p 
and A/c- 



2.3 More complicated density distributions 

Several complications may arise, including multiple pattern 
speeds, non-stationary pattern speeds, solar motion, etc. 
Here we discuss issues related to violations of the density 
assumption in Eqn. [l], deferring discussion of solar motion 
to 1 2.4. Possible violations of Eqn. |l| include: 



(i) Multiple pattern speeds. The MWG has a bar and spi- 
rals, perhaps also a small lopsidedness; almost certainly all 
have different pattern speeds. To understand how these will 
affect the measurement, we suppose that there are n non- 
axisymmetric features superposed, so that Eqn. [i] becomes 
p(x,y,t) = X/"-o wO 9 '^ 1 ~ ^«*)- Each of these asymme- 
tries gives rise to integrals V% and Id. Then, by the linear- 
ity of the continuity equation, we can write Ftp . f2jPj = 
J"\ /Q = /C. Dividing both sides by P = yj. Vi, we obtain 
RqQ — /C/P, where Q. is the average of the different pattern 
speeds, weighted by the Pi's. Thus the average is weighted 
by both the asymmetry and the density (since if [i\ = s/12 
then Pi — sVi). Therefore, small perturbations to the disk 
density do not contribute very significantly to fi. However, 
since the Pi's can be either positive or negative, Q, is not 
limited to values in the range min(f2i) and max(fli) (unless 
all the Pi's have the same sign). If distances are available, 
f{R) can be varied, to test for such cancellations; otherwise, 
strong changes in the slope of P as a function of lo may be 
used to test for cancellations. 

Note that winding spiral arms are a special case of the 
multiple pattern speed situation. In this case, Q p is a contin- 
uously varying function of distance from the Galactic center. 
This case is discussed further in § |6.l| . 

(ii) Radial oscillations. If such oscillations were present, 
Eqn. [i] would become p,(x,y,t) = p,{p + pie -1 "'-* - ' ', tp — 
0. v t), where we assume pi is small. The time-derivative term 
in the continuity equation therefore becomes 



dp 
dt 



-Q 



dp 

'di> 



dp 



OR 



(11) 



The ratio | Ai /^4o| is then given by ujpi/(2Q p Rd), where, for 
concreteness, we have assumed that the disk is exponential 
with scale-length Rd and the azimuthal density variation is 
of the form p(R)e 2t ^. For a simple estimate, we assume w ~ 
max(fc) ~ 2S1. Therefore we obtain A1/A0 ~ V c p\/£l P Rdp\ 
at p — 1 kpc, if the circular velocity is 220 km s _1 , Q p — 60 
km s _1 kpc -1 and Rd = 2.1 kpc (Bissantz et al. 2002), 
we obtain | Ai /j4q ~ 1.7pi kpc -1 . Thus radial oscillations 
with amplitude larger than about 100 pc will strongly in- 
terfere with the measurement of pattern speeds using the 
TW method. However, radial oscillations are expected to be 



strongly Landau damped, and it would be very surprising if 
the MWG exhibits large radial oscillations. 

(iii) Growing amplitudes. In this case, we can replace Eqn. 
[l] with p(x,y,t) = po(p) + pi(p,-ip — Q, p t)e at , where a > 
(q < 0) represents a growing (damped) feature. The term 
po(p) is included to make explicit the fact that the density 
of the feature grows at the expense of a background axisym- 
metric disc (for which P = K, = 0), to conserve the total 
mass. The time-derivative term is now replaced by 



dt ~ ilpe 



'^-+ap 1 e at = Aq + Al 

dip 



(12) 



It is possible to propagate the term in api through the cal- 
culation, giving, in the 2D case 



n p R V + a 



Rg(R)pi dl dR = K.. 



(13) 



Eqn. [l^ contains both f(R) and g{R), and therefore requires 
distance data for a solution. Lacking distances, we estimate 
the interference from growing features by the ratio |j4i/Ao |, 
which, for an m = 2 azimuthal density variation, is given 
by |a|/(2fi p ). Thus the growth term is important when the 
growth rate is of order the pattern speed (in which case it 
is safe to ignore damping). Since rapid growth is mostly 
restricted to the low amplitude, linear regime, while A-p 
can only be large enough for an accurate measurement once 
growth is nearly saturated, we expect that rapid growth is 
unlikely to significantly interfere with pattern speed mea- 
surements. 

The most important growing structures are spirals. If they 
are modal, they can grow even when tightly wound. The 
theoretically computed growth rates vary substantially; e.g. 
the fastest growing two arm modes in the flat-rotation, cut- 
out discs of Evans & Read (1998) have 0.03 <S a/2Q, p <, 0.29. 

(iv) Time-varying pattern speeds. In this case, Eqn. [I] gets 
replaced by p{x,y,i) — p(p,ip — Q p (t)t). Although this in- 
troduces a new term, proportional to t-gf-, it also changes 
the angle derivative to one in the bar rest frame. A simple 
transformation to a time il in which the moment of observa- 
tion is t' = recovers Eqn. |^. Thus a time-varying pattern 
speed does not disturb the TW method. 

We have considered several violations of the density 
condition of Eqn. |l| Of these, perhaps the most likely are 
growing spirals, but it is unclear how rapidly growing spi- 
rals are. Most likely, if spirals are of sufficient amplitude 
to give large signals, then their growth rate cannot be too 
large. This may be sorted with more detailed data on the 
MWG density distribution than we currently have, such as 
will be available with future astrometric missions. For the 
present work, we have assumed that there are no (strong) 
non-axisymmetric structures with rapidly growing ampli- 
tudes. 



2.4 Solar motion 

We derived Equation (^) for an inertial frame. Any radial 
velocities that can be measured will be heliocentric, and 
therefore the motion of the Sun needs to be taken into ac- 
count. The Sun's velocity can be decomposed into two parts: 
a motion of the local standard of rest, Vlsr, and the pe- 
culiar motion of the Sun relative to the LSR, V0 iPec . We 
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assume that the motion of the LSR consists of two parts, 
a tangential motion around the MWG centre with velocity 
vlsr = ulsrx + Vlsry, where x is in the direction towards 
the Galactic centre and y is the in-plane tangential direction. 
Then, if velocities in the heliocentric frame are labelled as vj, 
we can write v r j = (v/ + v ) f ; = v' r>i +VQ, pec fi +v L sr ■ 
?i = v ' r ,i + v ©,pec ■ fi + Vlsr sin U cos bi + u l .sr cos U cos h. 
Thus Equation (^) can be rewritten: 



AV 



flpRo — Vlsr 

J2i f( r i) sin h cos bi 
J2 t f{ r i) cos ^ cosbi 



-ULSR 



J2i /( r sin ^ C0S bi 



K, S 

^ -WLSR^- 



(14) 



This equation, which is the 3D form of the TW equation de- 
rived by Kuijken and Tremaine (1991), is used in the analysis 
of the OH/IR star sample below. It contains two quantities, 
Vlsr and ulsr which must be determined from other data. 
Here, Vlsr includes the circular velocity; Reid et al. (1999) 
measured Vlsr = 219 ± 20 km s _1 for Ro = 8 kpc, from 
the apparent proper motion of Sgr A*, while from a similar 
measurement Backer & Sramek measured Vlsr = 234 ± 7. 
The uncertainty in mlsr leads to a much larger uncertainty 
in AT/, largely because S depends on the cost's, which 
do not cancel across 1 = 0. The TW equation for exter- 
nal galaxies is similarly sensitive to the relative motion of 
the observer and the observed galaxy, but in that case re- 
flection symmetry is used to minimize this sensitivity. In 
the absence of distance information, no similar symmetry 
argument can be deployed in the MWG. The solar pecu- 
liar motion, V0 iPec , has been measured with HIPPARCOS 
data by Dehnen & Binney (1998). They found vo. pec = 
(U , Vb, Wo) = (10.00, 5.25, 7.17) ± (0.36, 0.62, 0.38) km s _1 . 
(These correspond to velocity components that are radially 
inwards, in the direction of Galactic rotation and vertically 
upwards, respectively.) In the absence of a radial motion of 
the LSR, we will then be measuring AV = £l p Ro — Vlsr, 
i.e. the negative of the solar velocity in the bar's rest frame. 

In the coordinate system we are working in, a star mov- 
ing away from the LSR has a positive v r ; inside the so- 
lar circle, such stars will primarily be at positive I. There- 
fore, if a positive AV is measured, this will correspond to a 
non-axisymmetric feature in the inner MWG with a pattern 
speed fl p > Qq, which must therefore be rotating in the 
same sense as the disc. 

The variance in AV due to the errors in the Sun's mea- 
sured peculiar velocity is given by: 



2 



2 (J2if( r i) COsliCOsbi ) 2 
a u ~T 2 + a v 

, 2 (E,/( r shift,) 2 
+°~w — 



-p2 



(15) 



(where we have assumed that errors in position measure- 
ments are insignificant). Note that, in this expression, the 
dominant term is the coefficient of a^ , since cos I is every- 
where positive in the range — ■§ < I < ^. 



3 TESTS 

We now test the TW method on some simple models. Our 
goal is not so much to test for specific structures expected 
in the MWG, since such an approach would say more about 
our models, which we would like to keep simple. Rather, 
our goal here is to learn how reliable the TW method is 
under various observational constraints, such as the number 
of discrete tracers, signal strength, restricted ranges of I, etc. 
Errors from sources specific to our application are discussed 
in §§ 

Models used to test the TW method need to have veloc- 
ity fields which preserve the assumed density distribution, 
i.e. the assumption of Eqn. (|l|) needs to remain valid. How- 
ever, because the continuity equation is purely kinematic, 
a model density distribution need not generate the poten- 
tial required by the kinematics. Here we test the method 
on models for a bar and a spiral, with application to the 
OH/IR stars in the ATCA/VLA OH 1612 MHz survey of 
Sevenster et al. (1997a,b & 2001) in mind. Since OH/IR 
stars are concentrated in the inner MWG, our tests assume 
that the non-axisymmetric structure is inside the Solar cir- 
cle. We assume throughout that measurements of I, b and v r 
for individual stars have only small errors, which we ignore, 
as is appropriate for this survey. 



3.1 Bar pattern speeds 

Tests of the method for a realistic bar can be obtained from 
TV-body simulations. We tested Equation for a selection 
of the bars in the fully self-consistent, TV-body simulations 
of Debattista & Sellwood (2000), who measured £l p from 
the time rate of change of the phase of the m = 2 Fourier 
component of the disc density. Here we present the results of 
several experiments using their "maximum disk" simulation 
at t = 650, which is shown in Fig. |l], when the bar was 
fast (TZ = 1.1 ± 0.2). At this time, a B = (2.6 ± 0.6)R d 
(where Ra is the disc scale-length), which is about twice 
the size of the MWG bar. The viewer therefore was placed 
at 4 < Ro/Rd < 5 (within the disc). The angle of the bar 
to the Sun-center line, ip, was varied in the expected range 
(20° < t[> < 30°). The angle lo, within which the integrals V 
and K, were evaluated, was varied in the range 30° < lo < 
60°. The simulation included 102, 000 disc particles; here we 
describe two experiments in which samples of 500 particles 
were selected. In the first experiment, particles were drawn 
using f(r) a Gaussian with f = 2.25R c i and oy = 0.7Rd- 
The second experiment used a cosine probability: /(r) = 



with R„ 



6R d for r < R max and /(r) 







^""^V S m „ 2 

otherwise. 

Fig. ^| plots the results, from which it is clear that the 
best results are obtained at large A-p and Ajc; since the 
Gaussian probability selects particles in a more limited ra- 
dial range, it often gives larger A-p and Aje and smaller 
errors. However, the cosine probability is likely to be more 
realistic for a flux-density-limited sample of tracers, since it 
preferentially selects nearby objects over distant ones. Fig. 
|2| show that, when A-p <; 0.15 and Ajc <; 0.15, the expected 
error is on average ~ 17%, and never worse than 40%. 

Variations in ip, lo and Ro mattered to the errors in AV 
only insofar as they changed A-p and Aje. For the Gaussian 
f(r) experiment, increasing Ro decreased A-p and A/c, with 
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Figure 1. Contours of projected disc density of the ./V-body 
barred system used in the tests. Contours are logarithmically 
spaced. This is the "maximum disk" simulation of Debattista &; 
Sellwood (2000) at t = 650. 

a corresponding increase in the fractional error. However, 
changes in if) (within our range) and lo did not substantially 
change A-p and Ajc, and the error is not overly sensitive to 
variations in these parameters. For the cosine probability 
distribution, A-p and Aje were insensitive to variations in all 
3 parameters. It is particularly reassuring that the errors do 
not depend on lo, provided that the bar is inside Iq. 



3.2 Spiral pattern speeds 

Stationary spiral models based on density wave theory 
are somewhat harder to construct. Since we have already 
demonstrated that the TW method works well for bars, 
which have larger non-circular motions than do spirals, a 
crude spiral model with purely circular orbits suffices for 
a spiral test. We therefore generated an exponential disc 
in the radial range 1.5i?d < r < 5Rd, consisting of two 
components, an axisymmetric component with particles on 
circular orbits at velocity 2.5 in some arbitrary units and 
a material (rather than wave) spiral part, which consisted 
of particles on circular orbits with constant Q p = 1. The 
spiral density varied tangentially as 1 + cos 2(9 — Oo), where 
8o — Q p t — ln(r) cot 7, to produce two-armed logarithmic spi- 
rals of pitch-angle 7, which we set to 20°. Both components 
were modelled by a Gaussian vertically, with equal scale- 
heights. The axisymmetric component accounted for | of 
the total mass, producing a spiral contrast consistent with 
that observed in external galaxies (Rix & Zaritsky 1995). 
The spiral component is shown in Fig. || 

Once we generated this system, we computed the inte- 
grals of Equation [] for a variety of orientations. In Fig. ^, we 
plot the errors at lo > 30° as before for the bar; in all cases 
the error is ^ 20% when Ap <; 0.15 and A/c <; 0.15. Since 
our rather crude spiral model included no random motions, 



Figure 3. The simple spiral system. The 10,000 particles shown 
are all from the spiral component; for clarity, the axisymmetric 
component has not been plotted. The observer, for the orientation 
of Fig. ^, is marked by a square, with the line segments indicating 
1 = 0°, ±30°, ±60° and ±90°. The dashed circle then has a 
radius r = 5Rd and the dotted circles have radius r±o> = 5±3Rd- 

we expect that more realistic errors will resemble those we 
found for the bar, provided signals are large. 

Our analysis of one orientation on a sample of 500 stars, 
selected with a Gaussian probability of r = 5Rd and ay = 
3Rd, is shown in Fig. |^; at lo = 45°, the error on AV is 
7.5% for this case. Note that A-p and A/c are both less than 
0.2 at lo = 45°, which is smaller than the signal seen in 
the OH/IR stars (§|E|). The results for this orientation bear a 
stri king resemblance to those of the OH/IR stars discussed 
in j fOf 

We ran a variety of additional tests of the TW method 
on our simple spiral system, with varying system orientation, 
radial probability function, number of particles selected, and 
reduced spiral mass fraction. These cases gave results con- 
sistent with those shown in Fig. [s| unsurprisingly, we found 
that the error decreases as the number of particles is in- 
creased. 



3.3 Synthesis of the tests 

We have shown that, for 500 stars, we expect to be able to 
measure AV with average errors of ~ 17%, and always less 
than 40%, for either a bar or a spiral, provided that signals 
are large (A-p and Ate 0.15). Furthermore, we found that 
limiting lo does not lead to large errors, provided that the 
full non-axisymmetric structure is included. 

In addition to the tests described above, we also ex- 
perimented with other parameters, with results generally in 
agreement with those presented here. The notable excep- 
tion was the case of discontinuous f(r)'s, which gave large 
errors. Such errors can be understood in terms of noisy sur- 
face terms in the integrals. Lacking distance information, 
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Figure 2. The fractional errors in the measurement of AV for a fast bar using 500 particles. The left panel is for the Gaussian probability 
and the right panel for the cosine probability. The symbols represent (absolute) errors as follows: circles 0.0 < <5(AV)/AV < 0.2, crosses 
0.2 < |<S(A10/AV| < 0.3, filled circles 0.3 < \5(AV)/AV\ < 0.4, triangles 0.4 < \8(AV)/AV\ < 0.6, diamonds 0.6 < \8(AV)/AV\ < 0.8 
and squares 0.8 < |tf(AV)/AV|. 



m - 
O - 




Figure 5. Errors in AV in the analyses of simple spiral mod- 
els. The sub-samples include different radial probability functions, 
orientations and 30° < l < 45°. Points have |5(AV)/AV| < 20% 
while crosses have |(5(AV)/AV| > 20%. Clearly, in all cases where 
A-p ^ 0.15 and A*; ^ 0.15, the expected error in AV is 20%. 

f(r) can only be a continuous function, corresponding to 
the detection probability. 



4 THE OH/IR STAR CATALOGUE 



4.1 OH/IR stars 

OH/IR stars are oxygen-rich, cool AGB stars in the super- 
wind phase (Renzini 1981) with typical M ~ 10" 5 M Q yr _1 . 
The wind outflow velocity, V c , is ~ 10 — 30 km s _1 , with 
V e larger in the more massive (younger) OH/IR stars. The 
mass-loss rate is large compared to the stellar mass (typi- 
cally, 1 to 6 Mq) and the superwind phase is therefore be- 
lieved to be very short-lived ~ 10 5-6 yr (Whitelock & Feast 
1993). Thus they are rare objects. On the other hand, they 
are old objects, with ages ~ 1 — 8 Gyr, and are therefore 
dynamically relatively relaxed. The dusty circumstellar en- 
velope which forms from the outflowing material absorbs the 
stellar radiation and re-emits it in the infrared, pumping OH 
masers (Elitzur et al. 1976) in the process. In optically-thick 
envelopes, the strongest line of this maser is at 1612.23 MHz 
which is conveniently insensitive to interstellar extinction. 
The Doppler shifted profile from the front and back of the 
thin OH shell permits easy identification of OH/IR stars, 
which makes possible an unbiased survey of these objects, 
tracing the kinematics of the inner MWG. A full review 
of the properties of OH/IR stars can be found in Habing 
(1996). 



4.2 The ATCA/VLA OH 1612 MHz survey 

We have used the data from the OH/IR star catalogue to 
search for pattern speeds in the inner MWG. The data were 
obtained in 3 surveys: the ATCA-bulge and ATCA-disc sur- 
veys (Sevenster et al. 1997a,b) at the Australia Telescope 
Compact Array (ATCA), covering -45.25° < I < 10.25°, 
\b\ < 3.25° and the VLA-disc survey (Sevenster et al. 2001) 
at the Very Large Array (VLA), covering 4.75° < I < 45.25°, 
\b\ < 3.25°. The ATCA region was covered uniformly with 
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Figure 4. The results of analysis with 500 particles drawn from the simple spiral model in the orientation of Fig. In the top row arc 
V and A-p, in the middle row K. and Afc, and in the bottom row is shown the number of particles, N part ; c i es , and AV obtained. The 
correct value of Ay is shown by the dotted line. 



1449 pointings. On the other hand, only 92% of the VLA 
region was surveyed, with 965 pointings out of 1053 giv- 
ing useful data. In all, the surveys produced 793 detections. 
The ATCA-bulge and VLA-disc surveys have a small over- 
lap region (4.75° < I < 10.25°). In this overlap region, 27 
sources were observed in both surveys, thus only 766 in- 
dependent sources were found. However, the overlap region 
also included detections in either only the VLA-disc or only 
the ATCA-bulge survey, which is, in part, a manifestation 
of the intrinsic variability of OH/IR stars. 



Details of the reduction technique for the 3 surveys 
are contained in Sevenster et al. (1997a,b) and Sevenster 
et al. (2001). Of the 793 detections, no outflow velocity 
was measured in 105 sources. Since such objects can be su- 
pergiants or even star-forming regions, rather than OH/IR 
stars, and don't have well-determined velocities, we have 
excluded them from our analysis. 



The use of two separate instruments for compiling the 
catalogue introduces a systematic difference in the complete- 
nesses of the two sides of the MWG. We corrected for the 
(empirically determined) completenesses, rewriting Equa- 



AV 



E p (Cp'EsepShi^ cos 6, 



(16) 



where is a sum over pointings, ^ s gp ^ s a sum over stars 
detected in pointing p and C p is the completeness of pointing 
p. The completenesses consist of two terms: C p = C Pia C Ptn . 
The first term, C Pl a, is a measure of the fraction of the sur- 
veyed area that is within some threshold radius. For this 
purpose, all VLA pointings, including the missing ones, were 
treated identically, and the Cp ja 's were calculated by Monte- 
Carlo integration. The threshold radii were set by requiring 
a constant primary beam response (PBR) threshold, equal 
in both telescopes, for each pointing centre. The PBR of the 
two instruments is shown in Fig. [| Table [l] lists several PBR 
levels, with the corresponding offset radii and values of C p ,a- 
The second term, C P: „, is a completeness due to the 
noise level, and varies from pointing to pointing. The Cp, n 's 
were obtained from the observed OH/IR star cumulative flux 
density distribution, N*, shown in Fig. ^, as the fraction of 
stars between the flux limits which are brighter than the 
detection threshold of each pointing. This is computed in 
the following way: 



tion (H): 



A^(max(/ m ^,7V 5 (Xp) < f < f max ) 

-^*(/miTi — f — fmax^) 



(17) 
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Table 1. The VLA and ATCA maximum offsets and the resulting values of C P: a for various PBR levels. 
PBR level VLA offset (') C p , a (VLA) ATCA offset (') C p , a (ATCA) 



0.5 
0.6 
0.7 



13.6 
11.7 
9.9 



0.656 
0.489 
0.347 



15.3 
13.3 
11.2 



0.812 
0.623 
0.450 



CC0.6 



r larcmm 




Figure 6. The PBR for the two instruments used in the OH/IR 
survey. The dashed (dot-dashed) line is the PBR for the ATCA 
(VLA). 



Here, fmin and f m ax are the flux density limits we use, a p 
is the noise in pointing p and N s is the multiplicative de- 
tection threshold (N s = 7, 4, 6 empirically determined for 
VLA-disc, ATCA-bulge, ATCA-disc respectively). We chose 
fmin — 0.16 Jy, since this is the highest of the lower absolute 
detection thresholds of the 3 surveys, and no upper limit on 
/. The noise level was more variable in the VLA-disc survey 
then in the ATCA surveys. In Fig. ^, we plot the resulting 
distribution of the C p ,„'s. 

The three surveys had different velocity coverages: the 
ATCA-bulge survey covered -280 km s _1 < v r < 300 
km s -1 , the ATCA-disc survey covered —295 km s _1 < v r < 
379 km s _1 , while the VLA-disc survey had full coverage in 
the ranee -200 km s" 1 < v r < 210 km s" 1 . Thus the three 
surveys had full common velocity coverage from -200 km s _1 
to 200 km s -1 . A considerable number of detections were ob- 
tained at \v r \ > 200 km s _1 , but these are mostly at |2| < 5°. 
One detection was obtained at 270 km s _1 in the ATCA-disc 
survey, but as this turned out to be a single-peaked object 
anyway, it was excluded from our analysis. 

Typical errors in positions are < 1" for the ATCA sur- 
veys and 2" for VLA-disc survey. These position errors in- 
troduce negligible errors in the value of AV, as can be ap- 
preciated from Equation ^| The velocity errors are of order 
1 km s" 1 for the ATCA surveys, and 2.5 km s 1 for the 
VLA-disc survey. The velocity errors produce larger errors 
in AV, although, as will be seen below, these are still quite 
small. 

The values of v r reported in Sevenster et al. (1997a,b 
and 2001) included a correction for an assumed solar pecu- 
liar motion of 19.7 km s" 1 towards RA = 18:07:50.3, Dec 
= +30:00:52 (J2000.0). For this study, we are adopting the 



Figure 7. The cumulative flux distribution of OH/IR stars in the 
ATCA/ VLA OH 1612 MHz survey. We used this flux distribution, 
defined from all sources with highest peak observed at PBR > 0.8, 
to compute the noise completenesses, Cp,n of the pointings. The 
lower limit is at 0.16 Jy while the upper limit is here set at the 
brightest observed OH/IR star (147.6 Jy). 



C n (ATCA) 



C n (VLA) 



Figure 8. The distribution of the C p>n 's. On the left are shown 
the ATCA pointings, and on the right are the VLA pointings. 
The effect of the higher noise level in the VLA-disc survey is 
easily apparent. 



solar peculiar motion found by Dehnen & Binney (1998); 
we have therefore had to transform from the Sevenster et al. 
frame to that of Dehnen & Binney. 

The mass, and therefore age and luminosity, of OH/IR 
stars correlate with V e . For example, Sevenster (2001) finds 
that OH/IR stars with V c = 13 (17.5) km s~ have a mass of 
roughly 1.7 (4)M and therefore an age of 1.8 (0.2) Gyr. We 
can anticipate, therefore, that we will need to make cuts on 
V c . However, the V e 's have been measured in units of half the 
channel-widths (2.27 km s _1 for the VLA-disc survey, and 
1.46 km s _1 for the ATCA surveys). To compensate for this 
effective binning, we have randomized the outflow velocities 
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V e (km/s) 

Figure 9. The dependence of [Z| max on V c , showing that only 
the nearest of the low V c OH/IR stars are visible in the survey. 
From this figure, we conservatively estimate that only at Ve > fO 
km s _1 are stars visible throughout the survey region < 45°). 
Each point is the average of the positive and negative I, while the 
error bars are defined as half the difference. 



within each bin. We did this in a number of ways with similar 
results once random re-samplings were performed; we have 
therefore relied on the simplest approach, namely to add to 
each V c a random value distributed uniformly between -1.2 
km s~ x and 1.2 km s _1 . 

The luminosity, L», of OH/IR stars increases with V a : 
simple models predict V a 4 tx L«, while more sophisticated 
models predict V c 3,3 oc L» (Habing et al. 1994). Thus a 
flux density-limited sample will probe to different distances 
for different V e . In Fig. M we plot the maximum longitude, 
^|max, of OH/IR stars versus V e , from which it can be seen 
that OH/IR stars of low V e are found only at small \l\. We 
interpret this as resulting from the ATCA/VLA OH 1612 
MHz survey being able to detect only the very nearest faint 
OH/IR stars of a distribution ending inside the solar circle. 
For the faint stars, the effect of small scale number density 
fluctuations is large, both because only a small portion of 
the disc is seen and because the number of such stars is low. 
We have therefore chosen to exclude all OH/IR stars with 
V c < 10 km s -1 . 

The TW method requires a sample of stars that form 
a relatively relaxed population; a convenient cutoff criterion 
would be at least 2.5 rotations at 6 kpc. This requires ages 
greater than 0.5 Gyr, or V c < 16 km s -1 . However, since we 
have had to randomize the outflow velocities we used only 
stars with V c < 15 km s _1 , to minimize contamination by 
very young stars. Thus our final range in outflow velocities 
is 10 km s _1 < V e < 15 km s , which is expected to cover 
a distance range of 4 to 10 kpc from the Sun. The lower 
limit removes roughly 10% of stars from the sample, while 
the upper limit removes a further 35%. 



5 RESULTS 

5.1 A pattern speed in the disc 

We computed V , /C, Ap and A/c for 5° < lo < 45° in steps 
of 1°. We set f min = 0.16 Jy, PBR = 0.5 and excluded 
OH/IR stars with |u r | > 280 km s~ . We tried both using 
the ATCA and the VLA in the overlap region. At each lo, the 
values reflect averages over 100 Monte-Carlo experiments 



containing slightly different samples. In these experiments, 
the V e 's were randomized in their bins (as described above) 
and Gaussian random errors of a v — 5 km s _1 added to the 
radial velocities. These results are presented in Fig. [l^. It 
is immediately clear that there is effectively no difference 
between using the ATCA or the VLA in the overlap region. 
Therefore, since the ATCA survey tends to be less noisy, 
we use it from here on. For lo < 25°, both V and K, are 
small and roughly constant, while A-p and A/c are noisy, 
suggesting they are strongly affected by Poisson noise. In 
the range 25° < lo < 33°, V , K., A-p and A/c all grow very 
rapidly. At lo ~ 33°, there is a sudden change in the slope 
of the number of included stars, which is associated with 
almost constant V , K, Ap and Aje at lo > 33°. In this region 
both Ap and A/c are <; 0.15, which our simple tests of 
suggest is large enough for moderate error levels. Averaging 
over lo > 35°, we obtain K./V = 263 ±44 km s _1 , where we 
have used one half of the difference between the maximum 
and minimum value for the error estimate since the values of 
K/V at different lo are not independent. Then, if ulsr = 0, 
Ro = 8 kpc and Vlsr = 220 km s _1 (which we assume for 
the remainder of this paper except where otherwise noted), 
Qp = 60 ±6 km s _1 kpc" 1 . 

Another way of presenting the data of Fig. [jij is shown 
in Fig. [ll|, where we plot the same data for lo — 45° but 
without the averaging of the various values of K. and V re- 
sulting from the Monte-Carlo reshufnings of V e . In the upper 
panel of this figure, we also plot lines of constant Q p (40, 50, 
60, 70 and 80 km s _1 kpc -1 ) assuming Ro — 8 kpc and 
Vlsr = 220 km s _1 for ulsr = and for ulsr = 5 km s . 
For the ulsr = case, the various re-samplings have aver- 
age Q p = 57 km s _1 kpc -1 , while in the case of ulsr = 5 
km s -1 , average Q p — 45 km s _1 kpc -1 , graphically demon- 
strating the sensitivity of the value of Q p obtained to a small 
radial LSR motion. The bottom panel also shows the effect of 
varying Vlsr and Ro, demonstrating that Q p increases with 
decreasing Ro and/or increasing Vlsr- However, we will not 
be concerned here with errors on Vlsr, Ro and mlsr; Vlsr 
enters only in the ratio Vlsr/^?o, which has been determined 
to 10% accuracy, Ro is determined to even higher accuracy, 
while the effect of ulsr is considered in §3.4. Note that the 



value of AV obtained this way is slightly different from that 
obtained by averaging over lo; since the values of AV at dif- 
ferent lo are not independent, we quote values of AV and 
Q p based on lo = 45° only from now on. 



5.2 Is it real? 

Before proceeding further, we present evidence that the sig- 
nal we have measured is real, and not a result of a systematic 
error in either our analysis or the OH/IR star surveys. Since 
/C and V measure differences, any systematic error in our 
calculated completenesses of the VLA or ATCA survey will 
be amplified into a signal. Several reasons lead us to believe 
that this has not happened. Firstly, the signals of Fig. |l(j 
first become evident at a very clear region in the disc, which 
is well correlated with an observed structure. In Fig. we 
present contour maps of the surface density of OH/IR stars 
in the survey. These clearly show the presence of very signif- 
icant over-densities localized at I ~ 25° and I ~ 32° . While, 
there are roughly equal number of stars at positive and neg- 
ative I, (despite the missing pointings, smaller PBHW and 
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Figure 10. The TW analysis for the OH/IR stars as functions of lo. In the top row are V and A-p, in the middle row K. and A/c, and 
in the bottom row is shown the number of stars, N s t ar s, and the resulting IC/V. In all panels, the solid line corresponds to using the 
ATCA survey in the overlap region, while the dashed line shows the result of using the VLA survey. In the bottom-right panel, circles 
(squares) are for the ATCA (VLA) survey, and the horizontal lines indicate the value of IC/V obtained by averaging results at Iq > 35°. 
The values of K. and the (ltr) error bars on K./V were obtained from 100 Monte-Carlo iterations as described in the text. 



greater noise in the VLA survey), for 25° < |/| < 33°, there 
are 13 OH/IR stars at negative I and 24 at positive I, which 
has less than 10% probability, even without taking the lower 
VLA detection probability into account. Thus, when all the 
Cp's are set to 1 (i.e. when no completeness corrections are 
included), the qualitative behavior of V, K,, A-p and Ak is 
unchanged, although, as is to be expected, signals are then 
smaller. 

Systematic errors in the completenesses are not likely 
to be the source of this signal. Suppose errors of the form 
Catca = Catca(1 + ^Catca) and/or Cy LA = CVla(1 + 
<5CVla) were present, and suppose further that the OH/IR 
stars sample only an axisymmetric part of the disc. It is 
then easy to show that <5CVla — 5C atca — V /V- which 
is approximately equal to —17/35 in our case. Since the 
noise errors are random, such a large systematic error would 
have to be in the values of the areal completenesses, which 
corresponds to roughly an overestimate (underestimate) by 
~ 20% in the radius at any given PBR for the ATCA (VLA). 
Such a large error is impossible since the PBR of both in- 
struments is well determined. A systematic effect might be 
introduced by the different channel widths of the VLA and 
the ATCA, but the similarity of the flux density distribu- 



tions between the VLA and the ATCA side suggests that, 
at most, this is a 10% effect. Furthermore, this effect acts to 
under-estimate the flux density of the sources detected by 
the VLA, so that f m in would remove from our sample more 
VLA stars than ATCA stars, which would make V negative, 
which is the opposite sign from what we found. 

Further evidence that the signal of Fig. [jl] is real comes 
from the dependence of the signal on b. We explored the ef- 
fects of limits in b, finding that most of the non-axisymmetric 
signal comes from very close to the MWG plane, with lit- 
tle signal at |6| > 1° (which divides our sample of OH/IR 
stars into two roughly equal parts). This is unsurprising for 
a disc population: at 6 kpc, 1° corresponds to roughly the 
disc scale-height. We found considerable fluctuations in the 
value of K./V at lo > 35° obtained from \b\ > 1°, reinforc- 
ing the impression that the settling of IC/V in the previous 
experiments is not caused by systematic errors, e.g. in the 
completenesses. 

All these facts make it likely, therefore, that the signal is 
intrinsic to the sample. Finally, it is perhaps worth drawing 
attention to the suggestive similarity between Fig. ^ and 
Fig. 0. 

Although the derived value of IC/V is nearly constant 
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Figure 11. The values of K, and V at Zo = 45° with bounds 10 < 
V c < 15 km s — 1 after reshuffling the V c 's. Each point represents 
a different reshuffling. Circles have A-p > 0.15 and A*; > 0.15, 
while triangles have either A-p or A/c or both less than 0.15. In the 
top panel, the dashed lines and dotted lines have Q p = 40, 50, 60, 
70 and 80 km s~ x kpc -1 (in order of increasing slope) assuming 
i?0 = 8 kpc and Vlsr = 220 km s — 1 . The dashed lines assume 
wlsr = 0, while ulsr = 5 km s — 1 for the dotted lines. The solid 
line (in both panesls) is the average Q p = 57 km s — 1 kpc -1 . In 
the bottom panel, we show the effect of changing Vlsr to 180 
km s _1 (dashed lines) and Rq to 7.5 kpc (dotted lines). 
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Figure 12. The surface density of OH/IR stars detected in the 
surveys. Adaptive smoothing has been used. The surface density 
has not been corrected for differences in the detection efficiencies 
between the ATCA (-45° < I < 10.25°) and the VLA (10.25° < 
I < 45°). Therefore the two sides of the plot (separated by the 
vertical dashed line) need to be considered separately. Note the 
very strong over-densities at 25° ^ I 5s 32°). The bottom panel 
shows the sample of OH/IR stars used in our analysis. 

in the range 35° < la < 45°, we cannot completely exclude 
that extending the survey to larger Zo would not lead to a 
different K/V. However, the density of OH/IR stars is de- 
creasing rapidly at Zo = 45° : an extrapolation of the binned 
surface density suggests that their surface density will van- 
ish at around Zo — 60°, as is also suggested by the dynam- 
ical models of Sevenster et al. (2000). Observationally, Le 
Squeren et al. (1992) found that the distribution of OH/IR 
stars is concentrated inside |Z| < 70°, with only roughly 15% 
of sources outside this range. Physically, this steep drop-off 
is the result of a metallicity gradient in the MWG disc: at 



the lower metallicities typical of the outer disc, AGB stars 
are more likely to evolve into carbon stars. Thus, even if 
there is still non-axisymmetric structure outside I — ±45°, 
it is not likely to change the derived value of K/V for the 
OH/IR population very significantly. 

5.3 Missing pointings and radial velocities 

The error introduced by the missing VLA pointings is sys- 
tematic in K and V, in that it always reduces the VLA (i.e. 
positive) contribution to these terms; however, the sign of 
the resulting change in K/V is not obvious. Of 88 missing 
pointings, 18 are in the overlap region. In general, therefore, 
we chose to use the ATCA pointings in the overlap region, 
to minimize the effect of missing pointings. To estimate the 
effect of the remaining missing pointings (70 of 88), we tried 
to reconstruct the missing pointings by reflection about the 
mid-plane where possible (49 pointings). Where this was not 
possible (21 pointings), we have simply averaged the value 
in the 4 nearest pointings, or the subset of them which were 
not themselves missing pointings. With these changes but 
all else as before, V, K, A-p and Ajc increase somewhat, 
but K/V changes to 239 ± 30 km s _ , which is not signif- 
icantly different from our previous values. The reason for 
this robustness to missing pointings can be understood by 
comparing the number of pointings and of OH/IR stars. 
Our sample contains 73 OH/IR stars detected in 840 good 
VLA pointings. With these probabilities, we estimate that 
70 x 0.09 ~ 6 OH/IR stars have been missed because of 
missing VLA pointings not in the ATCA overlap region. 

Due to observational constraints, we considered only 
OH/IR stars with \v r \ < 280 km s" 1 . We experimented with 
reducing the maximum |w r | further to test how our limit 
might have affected our value of K/V. Radial velocities this 
large are expected to occur only in the bulge region, and 
indeed, the biggest relative changes in K were at small Zo, 
but these propagated to small changes in K at large Zo. It 
is only at |v r | < 250 km s _1 that stars begin to leave the 
sample, and K/V decreases slowly for smaller maximum \v r \, 
reaching 229 ±40 km s" 1 at \v r \ < 200 km s" 1 , which is still 
smaller than la. We therefore conclude that it is unlikely 
that our modest cuts in \v r \ has resulted in a large systematic 
error in K/V. 

5.4 Radial motion of the LSR 

A non-zero radial velocity of the LSR, ulsr, cannot be ex- 
cluded. For our OH/IR star sample, S/V — 18. Thus even 
a quite small Ulsr will drastically alter the value of AV: if 
ulsr — 14 km s - , then AV changes by 100%. An accurate 
measurement of ulsr is therefore needed. 

Averaging the 197 OH/IR stars at -10° < Z < 10° (i.e. 
the ATCA-bulge survey), we estimated ulsr = ±2.7 ± 6.8 
km s _1 (the plus sign indicating that the motion is away 
from the Galactic centre). Similarly, Kuijken & Tremaine 
(1994) used a variety of tracers (OH/IR stars, globular clus- 
ters, high velocity stars, planetary nebulae, etc.) to obtain 
itLSR = — 1±9 km s _1 . Perhaps the best constraint on ulsr 
is that of Radhakrishnan & Sarma (1980), who showed that 
the main component in the HI absorption spectrum of Sgr A 
has a mean line-of-sight velocity 7/r = —0.23 ± 0.06 km s _1 . 
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ent from the results we presented above. However, we found 
that the average Q, p tends to systematically increase as the 
upper cutoff on V c is decreased. Thus, for example, we find 
that for 10 < V c < 14.5 km s" 1 (212 stars), the average 
tt p = 64 km s" 1 kpc" 1 while for 10 < V e < 14.0 km s" 1 
(181 stars), the average Q p = 69 km s _1 kpc" 1 . The same 
trend was seen when we increased the upper V e to 16 km s _1 , 
although in that case we are less certain that the popula- 
tion is relaxed. While all these results are consistent with 59 
km s" 1 kpc" 1 within 2a, the systematic behavior suggests 
that there is a systematic error, perhaps due to the popula- 
tion not being fully relaxed. We take as an estimate of this 
error 10 km s _1 kpc -1 . 



3> 



Figure 13. The distribution of values of K. and V in re-sampling 
experiments with 220, 230 and 240 OH/IR stars. The list of stars 
is first sorted by V e , then consecutive stars are sampled in the 
range 10 < V e < 15 km s — . The dashed lines indicate Q p = 40, 
50, 60, 70 and 80 km s" 1 kpc" 1 . 

Since this absorbing material is presumably outside the in- 
ner region of the MWG and within 6-7 kpc from the Sun, 
this means that a wide band of material is moving at a com- 
mon radial velocity. While it is possible to construct models 
in which dvn/dR is zero but with non-zero wlsr (e.g. Blitz 
& Spergel 1991), the most natural explanation is that mlsr 
is zero to this level; if ulsr = —0.23 km s - , this changes 
AV by -4 km s -1 . 

5.5 Sampling experiments and error estimates 

We now estimate the errors arising from sampling, which 
should be the largest source of errors for our small sample. 
We do this by re-sampling the data; however, since arbi- 
trary resamplings will increase the noise, we first sorted the 
OH/IR stars by V c , then took sub-samples of 220, 230 and 
240 consecutive stars. At each sub-sample size, we gener- 
ated 1000 different reshuffiings of the V c 's (since we always 
used 10 < V e < 15 km s , the number of combinations of 
OH/IR stars, and therefore the number of different values 
of K, versus V is smaller than 1000). The results are pre- 
sented in Fig. |l^, which shows clearly that there is quite a 
fair amount of scatter in the values of K, and V . The me- 
dian value of AV = 252 ± 41 km s - ; this corresponds to 
Sip = 59 ± 5 km s" 1 kpc" 1 . Therefore, as an estimate of the 
sampling noise, we adopt 5 km s -1 kpc -1 . 

The error estimate on AV from sampling is consistent 
with the results we found in and with the estimate based 
on variations in lo. We further test this error estimate by 
taking a sub-sample of PBR = 0.6, which gives about 180 
stars. Then the average AV = 175 ± 23 km s - , which is 
within 2a of the result for the full sample. Since this is an 
even smaller sample than our main sample, it should be even 
noisier than our full sample. We conclude that our sampling 
error error estimate is reasonable. (Higher PBR levels give 
even smaller samples and require larger areal corrections 
(see Table 0), which make them too susceptible to noise. We 
have therefore not attempted to use PBR levels above 0.6.) 

We have experimented with other sub-samples; in all 
cases we found that the results were not statistically differ- 



6 DISCUSSION AND CONCLUSIONS 

The results of our study of the OH/IR stars in the 
ATCA/VLA OH 1612 MHz survey give us a pattern speed 
of fl p = 59 ± 5 ± 10 km s -1 kpc" 1 (internal and estimated 
systematic errors respectively) for mlsr = 0, Ro — 8 kpc 
and Vlsr = 220 km s . If the rotation curve stays flat 
between the region of the non-axisymmetry and the Sun, 
then the corotation radius is 3.7 ± 0.31q g kpc, while the 
main non-axisymmetric feature in question appears to be 
centered at I ~ 30°, i.e. ~ 4.0 kpc if it is seen tangentially. 
If this is the full extent of this non-axisymmetric feature, 
then 1Z = 0.9 ± O.iIq'i, which is fast in the usual definition. 
However it is possible that the feature is not being viewed 
tangentially and extends somewhat further out, in which 
case 1Z would be smaller. 

Which feature in the MWG disc can be responsible for 
the signal we have measured? While the TW method has 
the great advantage of being model independent, the result 
is that, in the absence of distance information, it is hard to 
identify the feature responsible for the pattern speed mea- 
sured. For such a fast feature, the MWG bar would seem to 
be the natural explanation. However, the longitude at which 
the non-axisymmetric signal peaks is somewhat larger than 
what present models of the MWG would have (e.g. Gerhard 
2001, but see Hammersley et al. 2000, who have recently 
claimed evidence of a second, larger bar in the MWG). 

The small latitude of the main part of the signal implies 
a disc source, while the large longitude suggests the signal 
arises mostly from spirals, particularly at the tangent point 
of the Scutum spiral arm. The high value of Q p and the 
relatively small 1Z may hint at a coupling to the bar; perhaps 
it is even an inner ring, rather than a spiral arm. Such rings 
are often found in external barred galaxies (Buta 1995), and 
are elongated along the bar, so that they co-rotate with the 
bars which they contain. Such a ring has been postulated in 
the MWG (e.g. Sevenster & Kalnajs 2001). 

The pattern speed measured must, in fact, be a den- 
sity and asymmetry weighted average of all the non- 
axisymmetric structure in the survey region, which is known 
to include a bar, various spirals and perhaps a ring. Note 
that there is no evidence of substantial cancellation of V 
in Fig. [l^. For various sub-samples, we find some marginal 
evidence for multiple pattern speeds: when we divided our 
full sample into bright (flux density > 0.854 Jy) and faint 
(flux density < 0.854 Jy) sub-samples, we found Q p — 49 ±6 
km s" 1 kpc -1 and Q. p — 71 ± 14 km s _1 kpc -1 respectively, 
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together with some evidence of signals peaking at different 
lo. Since the bright (faint) sub-sample is, on average, ex- 
pected to be closer to (further from) the Sun and further 
from (closer to) the Galactic centre, the lower (higher) pat- 
tern speed is not unexpected. However, these sub-samples 
are even smaller than our already small full sample, and the 
differences are not statistically significant. Multiple pattern 
speeds can be better studied with larger samples of tracers. 



6.1 Future prospects 

In the future, it will be possible to improve on our measure- 
ment in a variety of ways. Our main limitation was the small 
sample of objects; larger samples will provide one means of 
improving on our measurement. At present, the ATCA/VLA 
OH 1612 MHz survey is the only large scale complete survey 
of a population satisfying the continuity equation. Kuijken 
& Tremaine (1991) applied the TW method to the Galactic 
HI, which covers the entire disc, but this is unlikely to sat- 
isfy continuity, particularly in the presence of a bar. How- 
ever, systematic surveys of other populations are now in 
progress (e.g. planetary nebulae Beaulieu et al. 2000; SiO 
masers Deguchi et al. 2001) or planned (Honma et al. 2000); 
these may be used for similar measurements. 

Accurate distance information will lead to an impor- 
tant refinement in the method. Future astrometric satellites, 
such as ESA's GAIA, will measure radial velocities together 
with distance. Not only will GAIA provide samples many 
orders larger than the ATCA/VLA OH 1612 MHz survey, 
but it will also be possible to drop the assumption of one 
pattern speed. Then, if the pattern speed is some function 
of Galacto-centric radius only, we can write 



J fl p (p)p(p,ip) cosf/jp 2 f(R)dipdp 



(18) 



where (p, ip) are Galacto-centric polar coordinates, R is the 
usual distance from the Sun, and for simplicity we have ig- 
nored both the third dimension and LSR motion. This is a 
Fredholm integral equation of the first kind for the unknown 
function fi p (p); techniques for the solution of such equations, 
including numerical methods, are well established. The main 
difficulty with a solution of this equation probably arises in 
regularising the solution at helio-centric radii R at which the 
integrals are small or vanish, as must happen for a spiral, 
but a detailed discussion of these issues is beyond the scope 
of this paper. Such data will make it possible to test theories 
of spiral structure in greater detail than has been possible 
to now. 



6.2 Conclusions 

We derived a 3D version of the Tremaine- Weinberg method 
for the Milky Way Galaxy. The method, as developed here, 
is based on the assumption that the density distribution 
can be expressed as in Eqn. jjj. Thus we have assumed one 
pattern speed only; if multiple pattern speeds are present, 
t hen an average pattern speed is measured, as discussed in 
Moreover, we have assumed that any rapidly growing 



structure present has low amplitude. We argued that this is 
likely to be the case since large amplitude features cannot 
sustain large growth for long. 



We tested the method on simple models of bars and 
spirals, selecting samples of 500 discrete tracers, to show 
that it is possible to measure pattern speeds with aver- 
age error of 17%, and always better than 40%, provided 
that asymmetry signals are sufficiently large, as described 
in We then extracted a sample of some 250 OH/IR stars 
from the ATCA/VLA OH 1612 MHz survey and applied the 
method to these stars. This gave K./V = 252 ± 41 km s _1 , 
from which we obtain fi p = 59 ± 5 km s _1 kpc~\ with 
a possible systematic error of perhaps 10 km s _1 kpc , 
if Vlsr = 220 km s -1 , Ro = 8 kpc and ulsr = (for 
other values of Vlsr, Ro or ulsr, flp can be obtained from 
flp.Ro — Vlsr = 252 ± 41 — 18ulsr). The sample is quite 
small and the detection is only significant at about the 90% 
level. Future larger samples will improve this situation; we 
sketch how distance data combined with the projected po- 
sition and line-of-sight velocity data used in our study will 
lead to a substantial improvement in our knowledge of the 
pattern speed(s) in the Milky Way Galaxy. 

This work was partly supported by grant # 20-56888.99 
from the Swiss National Science Foundation. We thank the 
anonymous referee for comments which helped improve this 
paper. 
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